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A central question is how the conformational changes of proteins affect their function and the inhibition of this function by drug 
molecules. Many enzymes change from an open to a closed conformation upon binding of substrate or inhibitor molecules. These 
conformational changes have been suggested to follow an induced-fit mechanism in which the molecules first bind in the open 
conformation in those cases where binding in the closed conformation appears to be sterically obstructed such as for the HIV-1 
protease. In this article, we present a general model for the catalysis and inhibition of enzymes with induced-fit binding mechanism. 
We derive general expressions that specify how the overall catalytic rate of the enzymes depends on the rates for binding, for the 
conformational changes, and for the chemical reaction. Based on these expressions, we analyze the effect of mutations that mainly 
shift the conformational equilibrium on catalysis and inhibition. If the overall catalytic rate is limited by product unbinding, we find 
that mutations that destabilize the closed conformation relative to the open conformation increase the catalytic rate in the presence 
of inhibitors by a factor cxp{A AGc / RT) where AAGc is the mutation-induced shift of the free-energy difference between the 
conformations. This increase in the catalytic rate due to changes in the conformational equilibrium is independent of the inhibitor 
molecule and, thus, may help to understand how non-active-site mutations can contribute to the multi-drug-resistance that has been 
observed for the HIV-1 protease. A comparison to experimental data for the non-active-site mutation L90M of the HIV-1 protease 
indicates that the mutation slightly destabilizes the closed conformation of the enzyme. 

Keywords: enzyme dynamics, induced fit, conformational selection, HIV-1 protease, non-active-site mutation, multi-drug 
resistance 
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1. Introduction 

The function of proteins often involves conformational changes 
during the binding or unbinding of ligand molecules |1, 2J- 
Central questions are how these conformational changes are 
coupled to the binding processes, and how they affect the func- 
tion of the proteins and the inhibition of this function by drug 
molecules. For some proteins, a conformational change has 
been proposed to occur predominantly after a binding or un- 
binding process OJjSl, apparently 'induced' by this process [7 1. 
For other proteins, a conformational change has been suggested 
to occur predominantly prior to a binding or unbinding process 
1(8 - 16 1, which has been termed 'conformational selection' since 
the ligand appears to select a conformation for binding or un- 
binding fTT-TPl. Binding via conformational selection implies 
induced-change unbinding, and vice versa, since the ordering of 
events is reversed in the binding and unbinding direction lfT6ll . 

In this article, we extend classical models of enzyme catal- 
ysis and inhibition [20] by including a conformational change 
during the binding and unbinding of substrate, product, or in- 
hibitor molecules. Our aim is to investigate how the confor- 
mational change affects the catalytic rates in the presence and 
absence of inhibitor molecules, and how non-active-site muta- 
tions that shift the conformational equilibrium alter these cat- 
alytic rates. We focus on enzymes with induced-change bind- 



ing mechanism since many enzymes close rather tightly over 
substrate or inhibitor molecules during binding. Binding via 
an induced-change mechanism, i.e. prior to the change from 
the 'open' to the 'closed' conformation of these enzymes, is 
required if the entry and exit of the ligand molecules are steri- 
cally obstructed in the closed conformation |3 |. In the reverse 
direction, the ligands then unbind via conformational selection 
because the change from the closed to open conformation has 
to occur prior to the unbinding process. Classical models of 
competitive inhibition, in which the inhibitor binds to the same 
site as the substrate, involve four states: an empty state E of 
the enzyme, and three states ES, EP, and EI in which the en- 
zyme is bound to a substrate molecule S, a product molecule 
P, or an inhibitor molecule I EOl . In our extended model for 
enzymes with induced-change binding mechanism, the enzyme 
can adopt two conformations, an 'open' conformation 1, and 
a 'closed' conformation 2. The model has seven states be- 
cause induced-change binding involves an open state EiS and 
a closed state £'26' with bound substrate, two such states EiP 
and E2P with bound product, and two states Eil and E2I with 
bound inhibitor, besides the empty state (see fig. 1). 

With our extended model, we derive general expressions 
for the catalytic rates of enzymes with induced-change binding 
mechanism. The catalytic rates of these enzymes depend on (i) 
the binding and unbinding rates of substrate, product, and in- 



Preprint submitted to BBA - Proteins and Proteomics 



February 19, 2013 



^ inhibitor binding ^ 

■ ^ b 



"12 



"21 



Ell 




substrate binding 



product unbinding 



E2S 



kr 



kf 



o 

a 
•< 

(A 




Figure 1: 7-state model for catalysis and inhibition of an enzyme with induced-fit binding mechanism. In this model, substrate 
molecules S and inhibitor molecules / first bind in conformation E\ of the enzyme. These binding events induce changes into 
conformation in which the substrate S is converted into the product P. In the case of the HIV-1 protease, the conformation E\ 
corresponds to the semi-open conformation, and the conformation E2 to the closed conformation. 



hibitor molecules in the open conformation 1, (ii) the forward 
and backward rates of the catalytic step, and (iii) the transi- 
tion rates between the two conformations in the bound states of 
the enzyme (see eqs. (j2]i to (jS}). Our general expressions for 
the catalytic rates lead to an effective four-state model with a 
single substrate-bound state, a single product-bound, and a sin- 
gle inhibitor-bound state (see fig. 2), but with effective on- and 
off-rates of substrate and product molecules that depend on the 
conformational transition rates (see eqs. (|6]l to ([8]l). 

The role of the conformational changes for catalysis and in- 
hibition can be revealed by non-active-site mutations that slightly 
shift the conformational equilibrium, but do not interfere di- 
rectly with binding and catalysis in the active site of the en- 
zymes. Several groups have suggested that such shifts in the 
conformational equilibrium might explain why non-active-site 
mutations can contribute to multi-drug resistance Il36ll42ll43l . 
i.e. to an increase of catalytic rates in the presence of different 
inhibitory drugs. Based on our general results for enzymes with 
induced-change binding mechanism, we investigate how these 
mutations affect catalysis and inhibition, and distinguish two 
cases. In case 1, the maximum catalytic rate of the enzyme is 
limited by the unbinding of the product. We find that the cat- 
alytic rate in the presence of inhibitors depends exponentially 
on the mutation-induced change AAGc of the free-energy dif- 
ference between the two conformations of the enzyme in this 
case (see eqs. (25i and (27i). Non-active-site mutations with 
AAGc > that slightly destabilize the closed conformation 
2 relative to the open conformation 1 of the enzyme lead to 
an increase in the catalytic rate, irrespective of the inhibitor. 
Such non-active-site mutations thus contribute to a multi -drug- 
resistance of the enzyme. In case 2, the maximum catalytic rate 
of the enzyme is limited by the forward rate of the catalytic 
step. In this case, mutation-induced changes of the conforma- 
tional equilibrium have no effect on the catalytic rate in the 



presence of inhibitors. A comparison with experimental data 
for the non-active-site mutation L90M of the HIV-1 protease 
indicates that this enzyme appears to follow case 1, which im- 
plies that non-active-site mutations that slightly destabilize the 
closed conformation contribute to multi-drug resistance. 

2. Catalysis and inhibition of an enzyme with induced-fit 
binding mechanism 

We consider an enzyme with two conformations E\ and 
E2 that binds its substrate S via an induced-change mechanism 
El ^ EiS ^ E2S. We assume that the catalytic step occurs 
in conformation 2 of the enzyme, and that the inhibitor / binds 
to the same site as the substrate ('competitive inhibition'). The 
catalytic cycle of the enzyme and the inhibition of this cycle 
then can be described by the 7-state model shown in fig. 1. 

The catalytic rate depends on the 14 forward and backward 
rates between the 7 states of the model. These rates are: 

(i) The binding rates s+[>5'], and and unbinding 
rates s^,p^, and i_ of substrate, product and inhibitor molecules 
in conformation 1 of the enzyme. Here, [S], [P], and [/] denote 
the concentrations of these molecules. 

(ii) The forward and reverse rates of the catalytic step k f and 
k'[' . 

(iii) The conformational transition rates S12, S21, P12, P21, ii2 
and 121 between the substrate-bound states EiSand E2S, product- 
bound states EiP and E2P, and inhibitor-bound states Eil 
and £2!- We assume that the bound states E2S, E2P, and 
E2I with conformation 2 of the enzyme are the ground-state 
conformations, while the conformations EiS, EiP, and Eil 
with conformation 1 of the enzyme are excited-state conforma- 
tions. This assumption is valid if experimental structures indi- 
cate that an enzyme adopts conformation 1 in its unbound state 
and conformation 2 in its bound states, since the experimental 
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Figure 2: Effective 4-state model for catalysis and inhibition of 
an enzyme with induced-fit binding mechanism. The effective 
binding and unbinding rates of substrate and product given in 
eqs. (|6]i to ([8| connect this model to the 7-state model of fig. 1. 

structures correspond to ground-state conformations. The as- 
sumption implies 



S21 < S12 , P21 < P12, and < ii2 



(1) 



i.e. the excitation rates S21, P21, and 121 are much smaller than 
the corresponding ground-state relaxation rates S12, P12, and 

One of our main results is that the catalytic rate can be writ- 
ten in the Michaehs-Menten form (see Appendix A) 



fcmax ['5'] 



(1 + [I]/Ki)K„-, + [S] 



with 



kfPoi 



and 



fc/ + kr + Pof 
{kf + 3of)pof + krSpf 
{kf + kr + Pof) Son 

S12S+ 



Sof ^ 



Pof 



S12 + S- 

S21S- 

S12 + S- 

P21P- 



(2) 



(3) 



(4) 



(5) 



(6) 



(7) 



(8) 



P12 +P- 

for negligible product concentration [P] or negligible product 
binding rate In deriving eqs. (j2|i to (jsjl from the ex- 

act solution, we have assumed that the excitation rate p2i in the 
product-bound state is much smaller than the relaxation rate S12 
in the substrate-bound state, besides eq. ([T]l (see Appendix A). 
This assumption is reasonable if the relaxation rates S12 and pi2 
in the substrate- and product-bound states are of similar magni- 
tude, since P21 P12 (seeeq. ([T]i) then implies also p2i ^ Si2- 



Our general results for the 7-state model of fig. 1 lead to 
an effective 4-state model of catalysis and inhibition, see fig. 2. 
The catalytic rate of this 4-state model is described by eqs. Q to 
(j4]| with Ki ~ «of/*on for negligible product concentration [P]. 
The effective 4-state model has the same structure as classical 
4-state models for competitive inhibition |20|, but includes two 
conformations for the enzyme. The rates Son and s^f in eqs. (j6| 
and (|7]i can be understood as the effective on- and off-rates of 
the substrate, i.e. as the effective forward and backward rates 
between the unbound state Ei and the bound ground state E2S. 
Similarly, pof in eq. (|8]l can be understood as the effective un- 
binding rate of the product from the bound ground state E2P. 
The expressions in eqs. (|6]l and (j8]l for the effective on- and 
off-rates result here from an analysis of steady-state catalytic 
rates, but are identical with previously derived expressions for 
the effective relaxation rates of two-step binding processes that 
involve a conformational change |fT6ll56l (see appendix B). 

3. Effect of non-active site mutations on catalysis 

In general, mutations can affect all rate constants of the 7- 
state model shown in fig. 1, in particular if they are located 
in the binding site, or active site, of the enzyme. Non-active- 
site mutations, in contrast, may mainly affect the conforma- 
tional equilibrium of the bound ground states and excited states. 
The conformational equilibrium constant S21/S12 of the two 
substrate-bound states EiS and E2S and the equilibrium con- 
stant P21/P12 of the product-bound states EiP and E2P de- 
pend on the free-energy energy differences between the states. 
These differences can be decomposed into the free-energy dif- 
ference AGc = G(i?2) — G{Ei) between the conformations 
and the differences in the binding free energy of substrate and 
product molecules in the two conformations of the enzyme. If 
a non-active-site mutation only affects the conformational free- 
energy difference but not the binding free-energy differences 
between the conformations, we have 



^21/^12 

S2l/s 



12 



P21M; 



exp(AAGc/i?T) 



exp{AAGc/RT) 



(9) 



(10) 



P21/P12 

where the prime indicates rates of the mutant. Here, AAGc = 
AGJ;^ — AGc is the mutation-induced change of the free-energy 
difference between the conformations. 

If the relaxation rates S12 and pi2 into the bound ground 
states are much larger than the corresponding unbinding rates 
s_ and p_ from the excited states, we obtain 
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(11) 
(12) 
(13) 



from eqs. (|6|l to ([8|l, which simplify to Son — s+, ~ S21S-/S12, 
and pof ~ P21P-/P12 in this case. We have assumed here that 
the non-active-site mutation does not affect the binding and 
unbinding rate constants s+, s_, and Ahernatively, the 
eqs. ( [TT| and ( [T3| ) follow from the eqs. (|6]l to ([8]l if we assume 
that the mutation-induced change of the conformational free- 
energy difference mainly affects the excitation rates, which im- 
plies S21/S21 — exp{AAGc / RT) and s'12 — si2 for the con- 
formational transition rates of the substrate-bound states, and 
P21/-P21 — exp{AAGc / RT) and p'12 — P12 for the transi- 
tion rates of the product-bound states lfT6ll . The non-active-site 
mutation thus changes the effective off-rates of substrate and 
product molecules by the same factorcxp(AAG'c/^7'), irre- 
spective of the binding energies of these molecules, while the 
effective on-rates remain the same. 

With eqs. ( [TT| to ( [T3] l, we can now determine the effect of 
a non-active-site mutation on the catalytic quantities k^^^ and 
K„i of the enzyme, which are given in eqs. (|3]l and We 
assume that the non-active-site mutation does not affect the for- 
ward and backward rates kf and kr of the catalytic step, and 
consider two cases: 

Case 1: kf much larger than Po/and s,,/ 

If the forward rate kf of the catalytic step is much larger 
than the effective off-rates of product and substrate, we obtain 
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kfp, 



'of 



K„ 



fc/Pof + fcrgpf 
{kf + fcr)Son 



(14) 



(15) 



from eqs. (|3]l and Q. The effect of the non-active-site mutation 
on fc„,ax and is then characterized by (see eqs. ( 1 1 1 to ( 13 1) 
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exp(AAGc/i?T) (17) 



for k', 



kf and ki 



kfPof + fcrSof 

kr- The non-active-site mutation thus 
changes fc^ax and Km by the same factor exp(AAG'c/-R7^)- 
The ratio k^^j^/Km of these quantities therefore is not affected 
by the mutation: 

(18) 
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K' 



Case 2: kf and kr much smaller than pof and Sgf 

If the rates kf and kr for the catalytic step are much smaller 
than the effective off-rates of substrate and product, we obtain 
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'.kf 

Sof 
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(19) 

(20) 



from eqs. ([3]) and (j4]). The effect of the non-active-site mutation 
then can be described by 
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^max 



and 



K' 
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K 



fof 
Sof 



exp(AAGc/i?T) 



(21) 



(22) 



The non-active-site mutation thus changes Km by the same fac- 
tor as in case 1, but does not affect the maximal catalytic rate 



From these equations, we obtain 



U 11(1 
max I m 

"-max /Km 



exp(-AAGc/i?T) 



(23) 



4. Effect of non-active-site mutations on inhibition 



To quantify the effect of a non-active-site mutation on the 
interplay of catalysis and inhibition, we focus now on substrate 
and inhibitor concentrations [/] and [S] with [/] ^ Ki and 
[I\Km/Ki ^ [S]. For such concentrations, a substantial frac- 
tion of the enzyme is bound to the inhibitor, and the catalytic 
rate of the enzyme is (see eq. (|2]i) 



kmaxKi [S] 
Krn[I] 



(24) 



and, thus, proportional to the inhibition constant Ki and in- 
versely proportional to the inhibitor concentration [/] . For [/] ^ 
Ki or [I]K„i/Ki ^ [S], in contrast, the catalytic rate of the en- 
zyme is independent of [/] and Ki (see eq. If the mutant 
also fulfills the conditions [/] > K^ and mK^/K- > [S] at 
the considered substrate and inhibitor concentrations, the ratio 
of the catalytic rates of mutant and wildtype is 



k' 



K'h' IK' 

-"^ i ""max I -^^in 
^max / -^ni 



(25) 



The expression on the right-hand side of this equation has been 
termed the 'vitality' of the mutant ETIl . The vitality depends 
on the catalytic quantities fcmax and Km and on the inhibition 
constant Ki of the wildtype and mutant enzyme. 

The inhibition constant Ki, defined in eq. (jSj, is the dis- 
sociation constant of the inhibitor and the enzyme. The con- 
stant depends on the free-energy difference G(i?2^) — G{Ei) 
between the bound ground state E2I and the unbound state 
El. This free-energy difference can be composed into the free- 
energy difference AGc = G{E2) — G{Ei) between the con- 
formations of the enzyme, and the binding energy AG/ = 
G{E2l) — G(i?2) in conformation 2 of the enzyme. The ef- 
fect of a mutation on Ki thus can be characterized by 



K 



= exp((AAGc + AAGi)/RT) 



(26) 



where K'i is the inhibition constant of the mutant, and AAGq 
and AAG/ are the mutation-induced changes of the confor- 
mational free-energy difference and the binding free energy in 
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conformation 2. Besides affecting AGc, some non-active-site 
mutations that contribute to inhibitor resistance of an enzyme 
have been suggested to have an indirect effect on the binding 
free energy AG/ by sUghtly changing the shape of the bind- 
ing pocket in a way that does not interfere with substrate bind- 
ing 122^-541. This suggestion concerns inhibitors with a bound 
shape that is not confined within the shape or 'envelope' of the 
bound substrate. 



Table 1: Catalytic efficiency of mutant L90M relative to wild- 
type HIV-1 protease (data from Maschera et al. B^ ) 



With eq. (26 1 and the results of the previous section, we 
can now determine the vitality of non-active site mutants. We 
consider again the two cases of the previous section: 

Case 1: kj much larger than pof and Sgf 

If the forward rate of the catalytic step is much larger than 
the effective off-rates of product and substrate for both wildtype 
and mutant, we obtain the vitality 



TC'h' IK' 



exp((AAGc + AAG/) /RT) (27) 



Ki A^max I K^ 

from the eqs. ( [T8] l and (26 1. The effect of the non-active-site 
mutation on the free-energy difference between the conforma- 
tion thus is reflected by the factor exp{/S./S.G c / RT) in the vi- 
tality, while the effect of the mutation on the binding-free en- 
ergy of the inhibitor in conformation 2 of the enzyme is re- 



flected by the factor exp(AAG//i?T). Eq. (27i indicates that 
the change of the conformational free-energy difference induced 
by a non-active-site mutation can have a significant effect on the 
catalytic rate of the mutant in the presence of the inhibitors. 

Case 2: kf and kr much smaller than pof and Sgf 

If the rates for the catalytic step are much smaller than the 
effective off -rates of substrate and product, we obtain the vital- 
ity 

K'h' IK' 



Ki /tjnax / 



exp(AAG//i?T) 



(28) 



from the eqs. ( |23) and ( |26) . In this case, the vitality is only 
affected by the change AAG/ in the binding-free energy of the 
inhibitor. 



5. Comparison with experimental data for the non-active 
site mutation L90M of the HIV-1 protease 

To illustrate our general results for the effect of non-active- 
site mutations on catalysis and inhibition of enzymes with induced- 
fit binding kinetics, we consider here the well-studied muta- 
tion L90M of the HIV-1 protease as an example. The HIV-1 
protease changes from a semi-open to a closed conformation 
during binding of substrate or inhibitor molecules. The con- 
formational change from the semi-open to the closed confor- 
mation has been suggested to occur after the binding process 
since binding in the closed conformation appears to be steri- 
cally prohibited 1251 l26i . The conformational dynamics of the 
HIV-1 protease has been studied in NMR experiments Il27lj29l . 
by pulsed EPR spectroscopy ll3Ql432ll . and in molecular dynam- 
ics simulations [33,-41 J . In the native state, the HIV-1 protease 



substrate 


k' 

max 

T- 

"^max 


Kn 
Km 


^max 1 -^m 
^max / 


SQNY-PIVQ 


1.3 


1.4 


0.9 


ARVL-AEAM 


1.6 


1.3 


1.3 


ATIM-MQRG 


1.4 


1.6 


0.9 


PGNF-LQSR 






0.4 


SFNF-PQIT 


1.3 


1.9 


0.8 


TLNF-PISP 


1.2 


1.1 


1.1 


AETF-YVDG 


1.4 


1.2 


1.2 


RKIL-FLDG 


1.2 


0.9 


1.5 


fluorogenic 


1.4 


2 


0.7 



is a homo-dimer of two identical subunits with 99 residues. The 
active site is located in a tunnel at the dimer interface, which 
opens and closes via the motion of two 'flaps', accompanied 
by motion in other regions of the protease. Inhibitors of the 
HIV-1 protease that bind to the active site play a central role 
in anti-AIDS therapies. However, the efficiency of these ther- 
apies is impaired by viral mutations that lead to inhibitor re- 
sistance. Mutations that are associated with resistance occur 
both in the active site of the HIV-1 protease and distal to the ac- 
tive site. While active-site mutations can directly interfere with 
catalysis and binding, the contribution of distal, non-active-site 
mutations to resistance is not fully clear. Several groups have 
suggested that non-active-site mutations of the HIV-1 protease 
such as L90M contribute to inhibitor resistance by shifting the 
equilibrium between the semi -open and closed conformation of 
the protease |36, 42 , 43 1. Alternative explanations for the mu- 
tation L90M include an indirect effect on active-site residues 
1 44-47 1, and a reduction of dimer stability |48| since the mu- 
tated residue 90 is located at the dimer interface. 

Maschera and co-workers |42| have systematically investi- 
gated the catalytic efficiency of the wildtype protease and the 
mutant L90M with 8 different peptide substrates and a fluo- 
rogenic substrate (see Table 1). The mutation L90M leads to a 
slight increase of fccat and K^. Averaged over all substrates, the 
ratios of the experimental values for the mutant and the wild- 
type are k'^Jk^^„ = 1.35 ±0.05 and K'JK^ = 1.43 ±0.14, 
while the average value for the ratio of k'^^JK'^ and k^^/Km 
is 0.98 ± 0.12. These results are consistent with eqs. ( [T6| ) to 
18 1, which indicate a change of /ccat and K^i 



by the same 

factor cxp{A AGc / RT) due to a non-active-site mutation that 
mainly affects the conformational equilibrium. From the aver- 
age values of A;,'„3x/fcniax and K'^^/K„i obtained from the data of 
Maschera et al., the change in the conformational free-energy 
difference AAGc due to the mutation L90M can be estimated 
as about 0.2 kcal/mol, which indicates a slight destabilization 
of the closed conformation of the HIV-1 protease relative to 
the semi-open conformation. However, other groups have re- 
ported a decrease of /cmax^ K„i, and kj„^/Km due to the mu- 
tation L90M for three different substrates P9l - |5T1 . In general, 
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Table 2: Ratio of dissociation constants K'jj and Kd for the 
mutant L90M and the wildtype of the HIV-1 protease (calcu- 
lated from data of Shuman et al. ll54l ) 



inhibitor 


K'n/KD 


Amprenavir 


2.73 ±0.67 


Indinavir 


2.85 ±0.71 


Nelfinavir 


1.99 ±0.70 


Ritonavir 


4.54 ± 1.68 


Saquinavir 


2.93 ±0.77 



the catalytic efficiency of the HIV-1 protease strongly depends 
on the experimental conditions, in particular on the pH and salt 
concentration 1421 l52l l53l . 

Shuman et al. f54] have measured the dissociation constants 
Kd and K'j^ of the wildtype and the mutant L90M for five dif- 
ferent inhibitors with biosensor methods (see Table 2). Within 
the experimental errors, the ratios K'jj/K]j of the dissociation 
constants for the five inhibitors agree with the average value 
3.0 ± 0.4 of these ratios. For a non-active-site mutation with 
a negligible change AAG/ of the binding free energy in the 
closed conformation of the protease, we expect a change of Kd 
by the inhibitor-independent factor cxp(AAGc/i?T) from eq. 
( |26l ) since Kjj is identical with the inhibition constant Ki. From 
the mean value K'j^ / Kd ~ 3, we obtain the estimate A AGc — 
0.65 kcal/mol for AAG/ ~ 0, which is somewhat larger than 
the estimate AAGc — 0.2 kcal/mol obtained from the changes 
of fcmax and K,ri reported by Maschera et al. (see above). 

Inhibition constants Ki and K[ for the wildtype and mutant 
L90M have been measured by several groups for a variety of 
inhibitors (see Table 3). The ratios K[/Ki differ for different 
inhibitors, which appears to indicate mutation-induced changes 
AAG/ in the binding free energy for the closed conforma- 
tion that depend on the inhibitor However, results by different 
groups for the same inhibitor also differ (see data in Table 3 for 
the inhibitors Indinavir and Saquinavir), possibly due to differ- 
ent experimental conditions such as pH or salt concentrations. 
Overall, the majority of the experimentally determined ratios 
K[/Ki for L90M in Table 3 are larger than 1, which may be 
interpreted to point towards a destabilization of the closed con- 
formation relative to the semi-open, i.e. towards positive values 
of AAGc. 

6. Conclusions 

We have considered here a general 7-state model for the 
catalysis and inhibition of enzymes with induced-fit binding 
mechanism. This model extends classical models for catalysis 
and inhibition [20 1 by including changes between two different 
conformations of the enzymes, and helps to understand how 
non-active-site mutations that shift the conformational equilib- 
rium can affect catalysis and inhibition. The induced-fit bind- 
ing mechanism considered here applies to enzymes such as the 
HIV-1 protease that close over ligands in the bound ground- 
state conformation in a way that sterically prevents ligand exit 



Table 3: Ratio of inhibition constants K[ and Ki for the mutant 
L90M and the wildtype of the HIV-1 protease 



inhibitor 


reference 


K'JK, 


TnHinavir 


1501 


4.1 ± 0.8 




0471 


1.3 ± 0.3 




Ii49l 


1 6 ± 03 


Nelfinavir 




9.0 ± 1.5 


Ritonavir 


|[50l 


3.2 ± 1.0 


Saquinavir 


Il42l 

ma 


21 ±4 
28 ±9 






3.0 ± 1.0 


L-735,524 


ma 


5.8 ± 1.7 


VX-478 


ma 


2.7 ±0.3 


AG- 1343 


ma 


3.5 ± 1.0 


ABT-538 


ma 


6.7 ±2.4 


Lopinavir 




1.0 ±0.6 


QF34 




0.9 ±0.4 


Darunavir (TMC-114) 


|47| 


0.14 ±0.04 



and entry in this conformation 13] |25] |26l . The unbinding of 
product and inhibitor molecules from these enzymes follows 
a conformational-selection mechanism since a conformational 
change from the closed ground-state conformation to an open 
or semi-open exited-state conformation is required prior to lig- 
and exit 1 16 1. Our approach can be generalized to enzymes with 
other binding mechanisms, e.g. to enzymes with a conformation- 
selection binding and induced-fit unbinding mechanism. 
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Appendix A. Catalytic rates 

We illustrate in this section how eqs. (|2]) to (|8]l for the cat- 
alytic rate of the 7-state model shown in fig. 1 can be derived 
from exact results. To avoid exceedingly long equations, we 
focus on the special case with inhibitor concentration [/] — 0. 
The derivation in the general case with [/] > is analogous to 
the derivation in this special case. 

For [/] = 0, the two states Eil and E2I have zero probabil- 
ity. The 7-state model thus reduces to a model with 5 states. The 
steady-state probabilities P{Ei), P{EiS), P{E2S), P{E2P), 
and P{EiP) of these 5 states can be calculated from the 5 equa- 
tions (see e.g. |55 1) 

P{Ei)+P{EiS)+P{E2S)+P{E2P)+P{EiP) = 1 (A.l) 

P{EiS)s^ + P(EiP)p^ - P{Ei)s+[S] = (A.2) 
PiEi)s+[S] + P{E2S)s2i - PiEiS){si2 + -s-) = (A.3) 
P{EiS)si2 + P{E2P)kr ~ P{E2S){s2i + kf) = (A.4) 
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P[E2S)kf+P{E^P)p^2 - P{E2P){kr+P2i) = (A.5) 



While eq. ( |A.l[ i ensures probability normalization, the eqs. (A. 2 1 
to ( |A.5| l are the flux balance conditions for the states Ei, EiS, 
E2S, and E2P in the steady state, in which the inward flux 
into a state equals the outward flux from this state. We have 
assumed that the product concentration [P] or product binding 



rate p+ [P] are negligible in eq. (A. 2 1. 

The catalytic flux k of the enzyme is defined as the steady- 
state flux along the cycle: 



k = P{E2S)kf - P{E2P)kr 



(A.6) 



From eqs. ( A. 1 1 to ( |A.6| l, we obtain the exact result for catalytic 
flux 

a[S] 



k 



b + c[S] 



(A.7) 



with 



a 
b 



s+Si2kfP2iP- (A.8) 

kfP2lP-{si2 +S-) + krS2lS-(j)i2 + P-) 

+ S21S-P21P- (A.9) 

fc/(P2lP- + {Pl2 +P2I +P-)si2)s+ 

+ {K{Pl2 +P-) +P2lP-){3l2 + S2l)s+(A.10) 



For P21 ^ P12 and S21 ^ S12 (see eq. ([T])), eq. (A. 10 1 can be 
simplified to 



fc/((P21 + Si2)p- +Pl2Sl2)s+ 
+ {kr{pi2 + P-) + P2lP-)si2S+ 



(A.ll) 



If we further assume p2i <C S12 (see text after eq. ([8]) for justi- 
fication), we obtain 



{{kf +fc,.)(pi2 +P-) + P2lP-)si2S+ 



(A. 12) 



The eqs. \A.%) and ( |A.12[ i now lead to the maximal catalytic 
flux 



kfP2iP^ 



C {kf + kr){pi2 + P-) + P21P- 



(A.13) 



This expression for the maximal catalytic flux is identical to 
eq. p I with pof given in eq. (j8]l. Similarly, the eqs. (A.9 1 and 
(A. 12 I lead to an expression for K„i ~ b/c that is identical to 
eq. ((4 1 with eqs. ^ to (|8]l. 

Appendix B. Effective rates for induced-fit binding 

We consider here the induced-fit binding process of an in- 
hibitor / to an enzyme that can adopt two conformations Ei 
and E2: 

El ^ Ell ^ E2I (B.l) 

i- 121 

In the forward direction of this process, the conformational change 
occurs after the binding event, apparently 'induced' by this event. 
In the reverse direction, the conformational change precedes the 



unbinding step. The unbinding of the inhibitor thus follows a 
conformational-selection mechanism. Here, i+ and i_ are the 
rate constants for the binding and unbinding of the inhibitor in 
conformation 1 of the enzyme, and ii2 and 121 are rates for 
the conformational exchange between the bound ground state 
E2I and the bound excited state Eil. We assume that the con- 
centration [/] of the inhibitor is much larger than the enzyme 
concentration, which implies pseudo-first-order kinetics with 
approximately constant inhibitor concentration and, thus, con- 
stant binding rate i+ [I]. 

The effective on- and off-rates of the two-step process ( |B.l| l 
can be determined from the dominant relaxation rate of this pro- 
cess 17611561 . Since the excitation rate 121 of the bound excited 
state Ell is much smaller than the relaxation rate ii2 into the 
bound ground state E2I, the effective off -rate from E2I to Ei 
is LI 6. .561 : 

^o^ ^ -. -— (B.2) 

«12 + «- 

The effective on-rate ion follows from the condition that the ef- 
fective equilibrium constant ion/iof of the two-step process has 
to be equal to the product (i+/i-)(«i2/*2i) of the equilibrium 
constants for the two substeps. With eq. ( |B.2[ ), this condition 
leads to 

ii2«+ 



ii2 + i- 



(B.3) 
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